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DISCOVERY OF GAMMA-RAY EMISSION FROM THE SHELL-TYPE SUPERNOVA REMNANT RCW 86 

WITH H.E.S.S. 
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ABSTRACT 

The shell- type supernova remnant (SNR) RCW 86, possibly associated with the historical supernova 
SN 185, with its relatively large size (about 40' in diameter) and the presence of non-thermal X-rays 
is a promising target for 7-ray observations. The high sensitivity, good angular resolution of a few arc 
qq I minutes and the large field of view of the High Energy Stereoscopic System (H.E.S.S.) make it ideally 

suited for the study of the 7-ray morphology of such extended sources. H.E.S.S. observations have 
fNj 1 indeed led to the discovery of the SNR RCW 86 in very high energy (VHE; E > 100 GeV) 7-rays. 

With 31 hours of observation time, the source is detected with a statistical significance of 8.5c and 
is significantly more extended than the H.E.S.S. point spread function. Morphological studies have 
been performed and show that the 7-ray flux does not correlate perfectly with the X-ray emission. 
The flux from the remnant is ~10% of the flux from the Crab nebula, with a similar photon index of 
about 2.5. Possible origins of the very high energy gamma-ray emission, via either Inverse Compton 
scattering by electrons or the decay of neutral pions produced by proton interactions, are discussed 
k/, \ on the basis of spectral features obtained both in the X-ray and 7-ray regimes. 

Subject headings: gamma-rays: observations - supernova remnants: general- supernova remnants: 
03 1 individual RCW 86 - H.E.S.S. 
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1. INTRODUCTION 

Shell-type supernova remnants (SNR) are widely be- 
lieved to be the prime candidates for accelerating cosmic 
ray protons and nuclei up to 10 15 eV. A promising way of 
proving the existence of high energy hadrons accelerated 
in SNR shells is the detection of very high energy (VHE; 
E > 100 GeV) 7-rays produced in nucleonic interactions 
with ambient matter. VHE 7-ray emission has been 
detected recently i n several shell-type SNRs, especially 
from C assiopeia A ([Aharonian et al.l (120011) ■ I Albert et al.l 
|2003)), RX .11713 .7-3946 (lAharon an et al.ll2007aft and 
RX J0852.0-4622 (|Aharonian et alj l2007bfh These two 
latest sources both show an extended morphology highly 
correlated with the structures seen in non-thermal X- 
rays. Although a hadronic origin is probable in the above 
cases (|e.g. Berezhko fc Volk II2006D. a lep tonic origin can 
not be ruled out (|e.g. Porter et al.ll2006l h 
Another young shell-type SNR is RCW 86 (also known as 
G315. 4-2.3 and MSH 14-63). It has a co mplete shell in . 
radio (iKesteyen fc Caswell I U98% . optical (jSmith 1119971) 
and X-ravs (|Pisarski et al.l 1984 ) , with a nearly circu- 
lar shape of 40' diameter. It received substantial atten- 
tion because of its possible associa tion with SN 185, the 
first historical Galactic supernova (|Clark fc Stephenson I 
Il977f l. However, conclusive evidence for t his connection 
is still missing: using optical observations. iRosado et al] 
(1996) found an apparent kinematic distance of 2.8 kpc 
and an age of ~ 10 000 years, whereas recent observations 
of the North-East part of the remnant with the Chan- 
dra and XMM-Newton satellites strengthen the case that 
the event recorded by the Chinese in 185 A D was a su- 
perno va and that RCW 86 is its remnant (tVink et alj 
2006). In this case, a distance to the SNR of ~ 1 kpc 
can be estimated for a standard Sedov evolution sce- 
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nario (|Bocchino et alj [2000) . The X-ray spectrum ob- 
tained with the Einstein satellite was first represente d 
by a two-tem perature pla sma model (|Winkler I Il978h . 
Then, RXTE (IPetre et alj|1999f) and ASCA ob servations 
(jBamba et al.l ([20001 ) . iBorkowski et all ([200 If )), with a 
wider spectral coverage, were used to resolve a non- 
thermal component in the X-ray spectrum which can be 
well described by a soft power-law with a photon index 
of ~ 3. The large-sca le de nsity gradien t acros s RCW 86 
(|Pisarski et all (|1984f ) and IClaas et al.l (|1989l )) possibly 
suggests that the northern part could be the shocked 
half of a very low-density wind bubble plus dense shell 
from the progenitor star, and to this extent it could well 
be similar to RX J1713.7-3946 and RX J0852.0-4622. 
In its southern part, RCW 86 contains an HII region. 
Apparently, the gas density in this HII region is rather 
high and spatially extended. Therefore, the SNR shock 
has swept over an extended high density region in the 
South, wit h consequent h i gh ra dio and thermal X-ray 
emissions (Bocchi no et al.l [2000). With a diameter of 
about 40', RCW 86 is one of the very few non-thermal X- 
ray emitting SNRs resolvable in VHE 7-rays. H.E.S.S., 
with its high sensitivity, its good angular resolution and 
its large field of view is ideally suited for morphology 
studies of such an extended object. 

Evidence for 7-ray emission from RCW 86 was found 
using the CANGARO O-II instrument, but no firm de- 
tection was claimed (jWatanabe et al.ll2003|) . Here, we 
present data on RCW 86 obtained with the full H.E.S.S. 
array between 2004 and 2007. 

2. H.E.S.S. OBSERVATIONS AND ANALYSIS METHODS 

H.E.S.S. is an array of four imaging Cherenkov 
telescopes located 1800 m ab ove sea level in the 
Khomas Highland in Namibia (jHinton 1 12004|) . Each 
telescop e has a tesselated mirror with an area of 
107 m 2 (|Bernl5hr et alj l2003f ) and is eq uipped with a 
came ra comprising 960 photomultipliers ([Vincent et al.l 
12001 covering a field of view of 5° in diameter. Due to 
the effective rejection of hadronic air showers with the 
stereoscopic imaging technique, the H.E.S.S. telescope 
system can detect point sources near zenith at flux levels 
of about 1% of the Crab nebula flux with a statistical 
significanc e of 5 a standard devi ation in 25 hours of ob- 
servation ([Ah aroman 

The shell-type SNR RCW 86 was observed between 2004 
and 2007 with the complete H.E.S.S. array. After stan- 
dard data quality selection and dead time correction, the 
resulting live time is 31 hours. The observations have 
been carried out at zenith angles ranging from 38° to 
53°. The data were taken using the wobble mode where 
the source is offset from the centre of the field of view, al- 
ternating between 28 minute runs in the positive and neg- 
ative declination or right ascension directions; the mean 
offset angle of the data set used in this analysis is 0.7°. 
The energy threshold of the system increases with zenith 
angle: for the observations presented here, the average 
threshold was 480 GeV. 

The data were calibrated using stand ard H.E.S.S. cal- 
ibratio n procedures, as discussed by lAharonian et alj 
(2004). The data were analyzed using a Hillas parame - 
ter based method as described in lAharonian et all (2005) 
with standard cuts, which include a minimum require- 
ment of 80 photo electrons in each camera image. Two 
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different backg round estimation p rocedures were used, 
as described in iBerge et al.l (|2007| ). For 2D image gen- 
eration and morphology studies, the ring background 
method was applied with a mean ring radius of 0.7°. As 
this method uses an energy averaged radial acceptance 
correction, the reflected-region background method was 
applied for spectral studies. In this second background 
subtraction procedure, OFF events were selected from 
the same field of view and in the same runs as the ON 
events by selecting the region symmetric to the ON re- 
gion with respect to the camera centre. As a cross-check, 
a second analysis chain, sharing only the raw data and 
using the "Combined Model" analysis (|de Naurois et al.l 

2005) , was also applied to the data. The two analysis 
methods yield consistent results. 

3. RESULTS 

A clear VHE 7-ray signal of 8.5<r standard deviation 
and 1546 ±183 excess 7-rays is detected from a circular 
region of 0.45° radius, centered on (aj20oo = 14 ,l 42 m 43 s , 
^72000 = — 62°28'48"). This integration region was cho- 
sen a priori on the basis of the X-ray data obtained with 
the ROSAT satellite and fully encompasses the SNR. Fig- 
ure Q] shows the VHE 7-ray excess map of the 1.6° x 1.6° 
region around RCW 86. The map has been smoothed 
with a Gaussian kernel with a er of 4.8' to suppress sta- 
tistical fluctuations on scales smaller than the H.E.S.S. 
point-spread function (PSF). The VHE 7-ray excess from 
RCW 86 is significantly extended beyond the PSF of the 
instrument, which is illustrated in the bottom left cor- 
ner of Figure [1] Contours of constant significance are 
superimposed in white at the 4, 5 and 6cr levels. An 
excess map has also been produced with the so-called 
"hard cuts" for better gamma hadron separation, which 
includes a stricter cut of 200 photo electrons on the im- 
age size compared to the "standard cuts" , and was found 
to be compatible with Figure [TJ The VHE emission 
shown in Figure [1] is suggestive of a shell- like morphol- 
ogy. To test this hypothesis, the brightness profile of a 
thick shell projected along the line of sight and folded 
with the H.E.S.S. point-spread function was fit to the 
unsmoothed excess map. As illustrated in Figure [IJ the 
best fit (x 2 / ndl = 233.1/220) is obtained with an outer 
radius of 24.43'±1.79^ tat , a width of 12.39' ±4.22^ tat and 
a centre of the shell at (a J2 ooo = 14' l 42 m 42.96 s ±14.1f tat , 
5j 20 oo = -62°26'41.6"±66.5£ at ). 

Figure [2] shows the radial profile of the VHE excess 
relative to the fitted centre. The fit of the radial 
profiles to the data points results in a chi-square per 
degree of freedom of x 2 / n df = 2.85/7 for a projected 
shell (determined by outer ring radius, ring width and 
absolute normalization) which is not significantly better 
than the fit of a projected uniformly-emitting sphere 
characterized by a ring radius and a normalization 
factor (x 2 /ndf = 5.43/8). Also visible in Figure [T] is 
an apparent deficit of 7-rays at the western part of 
the SNR. However, the azimuthal profile in Figure [3] is 
consistent with a constant and reveals that this dip is 
not significant (x 2 /ndf = 1.47/5). 

Figure [4] shows the 3-6 keV X-ray map of RCW 86 
obtained using six observations of the remn ant carried 
out b y the XMM-Newton satellite in 2006 (|Vink et all 

2006) and additional observations taken in 2007. The 




14 h 45 m 14 h 40 m 

Right Ascension (J2000) 



Fig. 1.—: H.E.S.S. 7-ray image of RCW 86. The 
map was smoothed with a Gaussian function with a 
^smooth = 4.8' to reduce the effect of statistical fluctua- 
tions. The linear color scale is in units of excess counts 
per arcmin 2 . White contours correspond to 4, 5, 6 a sig- 
nificance, obtained by counting gamma rays within 0.14° 
from each given location. The image inset in the bottom 
left corner indicates the size of a point source as seen 
by H.E.S.S., for an equivalent analysis, smoothing and 
zenith angles. The centre of the fitted shell, as discussed 
in the text, is marked by a black cross. The two solid 
green circles correspond to the inner and outer radii of 
this shell. 
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Fig. 2. — : Upper panel: H.E.S.S. radial profiles around 
the fitted centre of the SNR (a J20 oo = 14 fe 42 m 43 s , 5j 2 qoq 
= -62°26'42"). The solid line shows the result of a 
projected uniformly-emitting sphere smoothed with the 
H.E.S.S. point-spread function and fitted to the H.E.S.S. 
data. The dashed line corresponds to a projection of a 
thick and spherically symmetric shell. The dotted verti- 
cal line illustrates the extent of the region used for the 
azimuthal profile and for the spectral analysis. Lower 
panel: Radial profiles of the X-ray data (3-6 keV) from 
XMM-Newton. These data are background subtracted 
and smoothed to match the H.E.S.S. angular resolution. 
Additionally, the obtained excess profile was normalized. 

energy range was selected to avoid as much as possible 
contamination from line emission from the, in general, 
cool plasma (< 1 keV) of RCW 86. Potentially, the 3 
- 4 keV range could contain some contamination from 
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Fig. 3. — : Upper panel: H.E.S.S. azimuthal profile 
integrated over a region of 0.5° radius covering the SNR 
RCW 86. The azimuthal angle is calculated with re- 
spect to the fitted shell centre. 0° corresponds to the 
North part of the source and 90° to the East. The solid 
line shows the result of a fit of the data to a constant 
which yields a chi-square of 1.47 for 5 degrees of freedom. 
Lower panel: Azimuthal profiles of the X-ray data (3-6 
keV) from XMM-Newton. These data are background 
subtracted and smoothed to match the H.E.S.S. angular 
resolution. Additionally, the obtained excess profile was 
normalized. 

Ar and Ca lines, but no such line emiss ion is seen 
in th e available Chandra, XMM-Newton (|Vink et all 
120061 ) or Suzaku spectra (|Ueno et all l2007f T This 
map was obtained by first automatically cleaning the 
observations of > 3a excursions to the mean count rate, 
thus minimizing the background of the maps. Then, 
for each observation and for each of the three detectors 
(MOS1, MOS2, and PN), a background count rate in 
the 3-6 keV band was determined using a relatively 
empty region of the field of view. In the final stage, 
the background image was subtracted from the count 
rate image, and then corrected using the exposure maps 
obtained with the standard XMM-Newton SAS 7.1.0 
software (which includes vignetting correction), in order 
to obtain the background corrected map displayed in 
Figure |U An overall positional agreement with the 
H.E.S.S. contours derived from Figure [T] as well as a 
good compatibility between the outer radius of the 
7-ray emission (24.43' ± 1.79g tat ) and the extension 
of the X-ray emission can be observed. However, the 
emission peak apparent in the X-ray azimuthal profile 
is not visible in 7-rays (Figure [3]). Furthermore, the 
dip in surface brightness at the center of the remnant 
seems more pronounced in the X-ray radial profiles 
(Figure [5]). A more detailed comparison of the 7-ray 
and X-ray morphologies would require higher statistics 
than presently available, and hence will have to await 
future longer observations. 

For the spectral analysis, the source region (ON region) 
is defined by a circle of 0.5° radius centered on the best 
fit position of the shell, chosen to fully enclose the whole 
source. The radius of the extraction region is illustrated 
in Figure O The spectrum obtained (see Figure [5]) is 
well described by a power-law with a photon index of 
2.54 ± 0.12 sta t ± 0.20 sys and a flux normalisation at 1 
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Fig. 4. — : Excess contours of 7-ray emission (0.55, 
0.8, 1.05 7-rays per arcmin 2 Gaussian smoothed with 
"smooth = 4.8') superimposed on the background sub- 
tracted XMM-Newton EPIC (MOS/PN) 3-6 keV X-ray 
image of the remnant. 



TeV of (3.72 ± 0.50 sta t ± 0.8 sys ) x 10" 
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(x /ndf = 6.30/4). The integral flux in the energy range 
1 - 10 TeV is (2.34 ± 0.3 sta t ± 0.5 sys ) x 10" 12 cm^s" 1 , 
which corresponds to ~ 10% of the integrated flux of the 
Crab nebula in the same energy interval. No significant 
improvement is obtained by fitting a power-law with an 
exponential cut-off (x 2 /ndf = 2.96/3). If the fit range 
is restricted to energies below 10 TeV, a photon index 
of 2.41 ± 0.16 s tat ± 0.20 sys and a flux normalisation at 1 
TeV of (3.57 ± 0.5 stat ± 0.8 sys ) x lO-^cm^s^TeV- 1 
are determined (x 2 /ndf = 0.68/2), compatible with the 
fit of the SNR in the whole energy range. 
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Fig. 5.—: Differential energy spectrum of RCW 86, 
extracted from a circular region of 0.5° radius around 
the position (a J20 oo = 14 h 42 m 43 s , S J2 ooo = -62°26'42") 
adjusted to the H.E.S.S. data to enclose the whole source. 
The solid line shows the result of a pure power-law fit. 
The error bars denote la statistical errors; the upper 
limit (arrow) is estimated at the 2a level. The bottom 
panel shows the residuals to the power-law fit. Events 
with energies between 600 GeV and 60 TeV were used in 
the determination of the spectrum. 
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4. DISCUSSION 

There are two commonly invoked mechanisms for VHE 
7-ray production in young supernova remnants, inverse 
Compton (IC) scattering of high energy electrons off am- 
bient photons (leptonic scenario) and n° meson produc- 
tion in inelastic interactions of accelerated protons with 
ambient gas (hadronic scenario) . In such a hadronic sce- 
nario, a comparison between the expected thermal X-ray 
emission and the actually measured thermal emission has 
to await deeper observations in which one can better de- 
termined whether the TeV emission traces the denser, 
thermal X-ray emitting parts of the SNR, or is more 
closely correlated with the X-ray synchrotron emission 
from the remnant. 

The measured 7-ray spectrum from RCW 86, restricted 
to energies below 10 TeV, translates into an energy 
flux between 1 and 10 TeV of 8.6 x 10~ 12 erg cm" 2 s _1 . 
The X-ray spectrum of the whole remnant is mixed be- 
tween thermal and non-thermal emission. Assuming that 
the hard X-ray continuum originates from non- t herma l 
synchrotron emis sion as reported by iRho et al.l (|2002f l 
IVink et all (120061) and lUeno et all (|2007fl . the measure- 
ment made bv lPetre et al. I (|1999ft using RXTE data pro- 
vides an estimate of the total amount of non-thermal flux 
from RCW 86. They find that the spectrum is well fitted 
by a power-law of index ~ 3 and a flux normalization 
at 10 keV of 10~ 4 cm~ 2 s _1 kcV -1 , which extrapolated 
down to the 0.7 to 10 keV band leads to an integral flux 
of 2.1 x 10~ 10 crgcm~ 2 s _1 . In a leptonic scenario, as- 
suming that the 7-ray emission is entirely due to the IC 
process on cosmic microwave background photons, the 
ratio of the synchrotron power and IC power radiated is 
often used to constrain the magnetic field. For a power- 
law distribution of electron energies, Kj~ p , the general 
equation relating the synchrotron power (Ps) produced 
by electrons with Lorentz factors between ji ; x and 72, x 
and the IC power (Pic) radiated between -fijc and 727c 
can be expressed as follow: 

Ps Ub (jg ~ itx) m 
Pic U ph {l *y c - 7 3-£) 

where U p h and Ub are the energy density of the pho- 
ton field and the energy density of the magnetic field, 
respectively. It should be noted here that, for a fixed 
X-ray energy, ji t x and 72, x are inversely proportional 
to the square root of the magnetic field. If X-rays and 
7-rays probe the same region of the electron spectrum, 
one finds the standard relation between the synchrotron 
and IC power -p^ = . Assuming that the target pho- 
ton field is the cosmic microwave background, a magnetic 
field of 30 /iG can be estimated using Equation [T] and the 
synchrotron photon index of ~ 3, independent of the dis- 
tance and a ge of the SNR. Th is estimate is compatible 
with that of IVink et alj (|2006l ) based on thin filaments 
resolved by Chandra (assuming a distance of 2.5 kpc) 
in which the authors also deduce a high speed of the 
blast wave (~ 2700 km s^ 1 ); their estimated value would 
increase to ~ 50 fiG for a distance of 1 kpc. However, 
it is still a factor of 2 lower than the m aximum field 
strength determined bv lVolkeF al. (2005) using a lower 
shock velocity of 800 km s -1 as sugg ested by optical dat a 
in the Southern region of the SNR (|Rosado et al.lll996l ). 
The difference between the field amplification estimated 



bv IVink et all (pOOl and that of iVolk et all ((20051) lies 
in the fact that Volk et al. obtained a higher result when 
they de-projected the measured filament width, as for an 
ideal spherical shock, whereas Vink et al. did not. With- 
out de-projection the two results remarkably agree, even 
though they were obtained for the southern side and the 
northern side, respectiv ely. A discu s sion o f de-projection 
for R CW 86 is given inlVolk et all (|2005h . With similar 
data, iBamba et al.1 (|2005l ) deduced a significantly lower 
magnetic field strength of ~ 4 — 12 /iG. However, their 
analysis is based on rather different assumptions on the 
nature of filament formation. 

In a hadronic scenario, one can estimate the total energy 
in accelerated protons W p in the range 10 — 100 TeV 
required to produce the 7-ray luminosity L 7 observed 
by H.E.S.S. using the relation Wp(10 - 100 TeV) « 

r 7 xL 7 (l-10TeV), in which r 7 « 4.9xl0 15 (r^y 1 s 
is the characteristic cool ing time of protons through the 
7T° production channel (jKelner et al.ll2006T ). The total 
energy injected in protons is calculated by extrapolating 
the proton spectrum down to 1 GeV. Because of this ex- 
trapolation over 4 decades in energy, the uncertainty of 
the estimate can be as large as a factor of 10. Assum- 
ing that the relatively steep slope of the proton spec- 
trum (as inferred from the observed 7-ray spectrum) is 
the result of an energy cut-off (somewhere around sev- 
eral tens of TeV in proton energy), and that at lower 
energies the proton spectrum has a E~ 2 type spectrum 
representative of those predicted by the diffusive shock 
acceleration theory, the total energy budget in all protons 
for the distance of 2.5 kpc and t he ambient gas densit y 
between 0.3 cm -3 and 0.7 cm - 3 (jBocchino et all [20001. 
would be (2 — 4) x 10 50 erg. This estimate is in rea- 
sonable agreement with theoretical expectations that a 
significant fraction of the explosion energy of 10 51 erg is 
released in relativistic protons. On the other hand, if 
the power-law spectrum of protons continues to GeV en- 
ergies with the spectral index T = 2.4 (i.e. similar to 
the gamma-ray spectrum below 10 TeV), the total bud- 
get in protons would exceed a few times 10 51 erg for a 
distance of 2.5 kpc. This would exclude the hadronic 
origin of TeV 7-rays, unless the SNR is nearby (~ 1 
kpc), o r the 7-rays ar e prod uced in very den s e re gions. 
Indeed. IPisarski et all (|1984h and lClaas et all (| 19891 ) re- 
ported that there is a large density contrast across the 
remnant, e.g. in the South, where the density could be 
as high as 10cm~ 3 ; with such a dense medium, a larger 
distance for the remnant could still be compatible with 
the observed 7-ray flux. 

5. CONCLUSIONS 

H.E.S.S. observations have led to the discovery of the 
shell-type SNR RCW 86 in VHE 7-rays. The 7-ray signal 
is significantly more extended than the H.E.S.S. point- 
spread function. The possibility of a shell-like morphol- 
ogy was addressed, but cannot be settled on the basis 
of the limited statistics available at the moment. The 
flux from the remnant is ^10% of that from the Crab 
nebula, with a photon index of about 2.5. The question 
of the nature of the particles producing the 7-ray signal 
observed by H.E.S.S. is also discussed. 
In a leptonic scenario, assuming that the 7-ray emission 
is entirely due to the IC process on cosmic microwave 
background photons and that the synchrotron and IC 
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photons are produced by the same electrons, the ratio of 
the 7-ray energy flux and the X-ray flux determines the 
magnetic field to be close to 30 fiG. 
In the hadronic scenario, the lack of information about 
the low-energy 7-ray spectrum results in large uncertain- 
ties on the total energy budget in protons. If below sev- 
eral tens of TeV, the proton spectrum has a E~ 2 type 
spectrum, the total energy in protons would be in rea- 
sonable agreement with theoretical expectations. On the 
other hand, if we assume that the proton spectrum con- 
tinues down to GeV energies with the observed spectral 
index T = 2.4, energetics would rule out a hadronic ori- 
gin for the TeV 7-rays unless the SNR is nearby, or if the 
7-rays are produced in a very dense medium as reported 
in the southern part of the remnant. 
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